Millimeter wave-based patient setup verification and motion tracking during radiotherapy.

BACKGROUND: Position verification and motion monitoring are critical for safe and precise radiotherapy (RT). Existing approaches to these tasks based on visible light or x-ray are suboptimal either because they cannot penetrate obstructions to the patient's skin or introduce additional radiation exposure. The low-cost mmWave radar is an ideal solution for these tasks as it can monitor patient position and motion continuously throughout the treatment delivery. PURPOSE: To develop and validate frequency-modulated continuous wave (FMCW) mmWave radars for position verification and motion tracking during RT delivery. METHODS: A 77 GHz FMCW mmWave module was used in this study. Chirp Z Transform-based (CZT) algorithm was developed to process the intermediate frequency (IF) signals. Absolute distances to flat Solid Water slabs and human shape phantoms were measured. The accuracy of absolute distance and relative displacement were evaluated. RESULTS: Without obstruction, mmWave based on the CZT algorithm was able to detect absolute distance within 1 mm for a Solid Water slab that simulated the reflectivity of the human body. Through obstructive materials, the mmWave device was able to detect absolute distance within 5 mm in the worst case and within 3.5 mm in most cases. The CZT algorithm significantly improved the accuracy of absolute distance measurement compared with Fast Fourier Transform (FFT) algorithm and was able to achieve submillimeter displacement accuracy with and without obstructions. The surface-to-skin distance (SSD) measurement accuracy was within 8 mm in the anterior of the phantom. CONCLUSIONS: With the CZT signal processing algorithm, the mmWave radar is able to measure the absolute distance to a flat surface within 1 mm. But the absolute distance measurement to a human shape phantom is as large as 8 mm at some angles. Further improvement is necessary to improve the accuracy of SSD measurement to uneven surfaces by the mmWave radar.

Script-based implementation of automatic grid placement for lattice stereotactic body radiation therapy.

Background and purpose: Spatially fractionated radiation therapy (SFRT) has demonstrated promising clinical response in treating large tumors with heterogeneous dose distributions. Lattice stereotactic body radiation therapy (SBRT) is an SFRT technique that leverages inverse optimization to precisely localize regions of high and lose dose within disease. The aim of this study was to evaluate an automated heuristic approach to sphere placement in lattice SBRT treatment planning. Materials and methods: A script-based algorithm for sphere placement in lattice SBRT based on rules described by protocol was implemented within a treatment planning system. The script was applied to 22 treated cases and sphere distributions were compared with manually placed spheres in terms of number of spheres, number of protocol violations, and time required to place spheres. All cases were re-planned using script-generated spheres and plan quality was compared with clinical plans. Results: The mean number of spheres placed excluding those that violate rules was greater using the script (13.8) than that obtained by either dosimetrist (10.8 and 12.0, p < 0.001 and p = 0.003) or physicist (12.7, p = 0.061). The mean time required to generate spheres was significantly less using the script (2.5 min) compared to manual placement by dosimetrists (25.0 and 29.9 min) and physicist (19.3 min). Plan quality indices were similar in all cases with no significant differences, and OAR constraints remained met on all plans except two. Conclusion: A script placed spheres for lattice SBRT according to institutional protocol rules. The script-produced placement was superior to that of manually-specified spheres, as characterized by sphere number and rule violations.

Inter-fractional portability of deep learning models for lung target tracking on cine imaging acquired in MRI-guided radiotherapy.

MRI-guided radiotherapy systems enable beam gating by tracking the target on planar, two-dimensional cine images acquired during treatment. This study aims to evaluate how deep-learning (DL) models for target tracking that are trained on data from one fraction can be translated to subsequent fractions. Cine images were acquired for six patients treated on an MRI-guided radiotherapy platform (MRIdian, Viewray Inc.) with an onboard 0.35 T MRI scanner. Three DL models (U-net, attention U-net and nested U-net) for target tracking were trained using two training strategies: (1) uniform training using data obtained only from the first fraction with testing performed on data from subsequent fractions and (2) adaptive training in which training was updated each fraction by adding 20 samples from the current fraction with testing performed on the remaining images from that fraction. Tracking performance was compared between algorithms, models and training strategies by evaluating the Dice similarity coefficient (DSC) and 95% Hausdorff Distance (HD95) between automatically generated and manually specified contours. The mean DSC for all six patients in comparing manual contours and contours generated by the onboard algorithm (OBT) were 0.68 ± 0.16. Compared to OBT, the DSC values improved 17.0 - 19.3% for the three DL models with uniform training, and 24.7 - 25.7% for the models based on adaptive training. The HD95 values improved 50.6 - 54.5% for the models based on adaptive training. DL-based techniques achieved better tracking performance than the onboard, registration-based tracking approach. DL-based tracking performance improved when implementing an adaptive strategy that augments training data fraction-by-fraction.

A Dose Accumulation Assessment of Alignment Errors During Spatially Fractionated Radiation Therapy.

PURPOSE: Spatially fractionated radiation therapy (SFRT) techniques produce high-dose peaks and low-dose valleys within a tumor. Lattice stereotactic body radiation therapy (SBRT) is a form a SFRT delivered across 5 fractions. Because of the high spatial dose gradients associated with SFRT, it is critical for fractionated SFRT patients to be aligned correctly for treatment. Here we investigate the dosimetric effect of daily alignment uncertainty through a dose accumulation study. METHODS AND MATERIALS: Dose accumulation was retrospectively performed for 10 patients enrolled on a phase 1 trial. Lattice stereotactic body radiation therapy was completed in 5 fractions with 20 Gy prescribed to the entire tumor and a simultaneous integrated boost of 66.7 Gy prescribed to a set of regularly spaced high-dose spheres. Daily alignment error was quantified through manually selected landmarks in both the planning computed tomography scan and daily cone beam computed tomography. The dosimetric effect of alignment errors was quantified by translating the isocenter in the treatment planning system by the daily average alignment error. Large errors were simulated by translating isocenter 5 and 10 mm for 1 and 2 fractions, independently assessing errors in the superior-inferior and axial directions. The reduction of dose gradients was quantified using the dose ratio (DR) of the mean dose in the high-dose and low-dose spheres. RESULTS: The average alignment error was 1.8 mm across the patient population resulting in minor smoothing of the high- and low-dose distributions in the dose accumulation. Quantitatively, the DR decreased from 3.42 to 3.32 (P = .093) in the dose accumulation study. The simulated worst case was an inferior-superior shift of 10 mm for 2 fractions where the average DR decreased to 2.72 (P = .0001). CONCLUSIONS: The dose accumulation study revealed on average DR only decreased from 3.42 to 3.32. However, setup errors >5 mm resulted in larger dosimetric degradation, reflecting a larger effect for individual high-dose spheres within regions exhibiting larger displacements.

Simulation study of a novel small animal FLASH irradiator (SAFI) with integrated inverse-geometry CT based on circularly distributed kV X-ray sources.

Ultra-high dose rate (UHDR) radiotherapy (RT) or FLASH-RT can potentially reduce normal tissue toxicity. A small animal irradiator that can deliver FLASH-RT treatments similar to clinical RT treatments is needed for pre-clinical studies of FLASH-RT. We designed and simulated a novel small animal FLASH irradiator (SAFI) based on distributed x-ray source technology. The SAFI system comprises a distributed x-ray source with 51 focal spots equally distributed on a 20 cm diameter ring, which are used for both FLASH-RT and onboard micro-CT imaging. Monte Carlo simulation was performed to estimate the dosimetric characteristics of the SAFI treatment beams. The maximum dose rate, which is limited by the power density of the tungsten target, was estimated based on finite-element analysis (FEA). The maximum DC electron beam current density is 2.6 mA/mm2, limited by the tungsten target's linear focal spot power density. At 160 kVp, 51 focal spots, each with a dimension of [Formula: see text] mm2 and 10° anode angle, can produce up to 120 Gy/s maximum DC irradiation at the center of a cylindrical water phantom. We further demonstrate forward and inverse FLASH-RT planning, as well as inverse-geometry micro-CT with circular source array imaging via numerical simulations.

A stochastic control approach to intrafraction motion management in intensity-modulated radiotherapy.

Objective.The goal of this research is to demonstrate proof-of-principle for managing intrafraction motion via feedback control of delivered dose to achieve dosimetry comparable to respiratory gating without compromising delivery efficiency.Approach. We develop a stochastic control approach for step-and-shoot intensity-modulated radiotherapy (IMRT) in which the cumulative delivered dose and future trajectory of intrafraction motion are dynamically estimated by combining pre-treatment four-dimensional computed tomography imaging and intrafraction respiratory-motion surrogates. The IMRT plan is then re-optimized in real time to ensure delivery of the planned dose in the presence of free-breathing motion. We compare the performance of the proposed approach against traditional motion-management techniques, namely, respiratory gating and internal target volume (ITV) planning, using the four-dimensional extended cardiac-torso computational phantom.Main results.We simulate the delivery of treatment plans for a lung tumor in the presence of variable breathing amplitude, tumor size, and location. Results show that the proposed method reduces irradiated tissue volume compared to ITV treatment. Additionally, it significantly reduces treatment time compared to traditional respiratory-gated treatment, without compromising the dosimetric quality.Significance.Respiratory gating is a common technique to manage intrafraction motion. While gating supports reduced treatment volumes, it also prolongs the treatment delivery time. The proposed stochastic control approach can help improve the delivery efficiency of respiratory gating without compromising the dose quality.

Phantom-based Quality Assurance of a Clinical Dose Accumulation Technique Used in an Online Adaptive Radiation Therapy Platform.

Purpose: This study aimed to develop a routine quality assurance method for a dose accumulation technique provided by a radiation therapy platform for online treatment adaptation. Methods and Materials: Two commonly used phantoms were selected for the dose accumulation QA: Electron density and anthropomorphic pelvis. On a computed tomography (CT) scan of the electron density phantom, 1 target (gross tumor volume [GTV]; insert at 6 o'clock), a subvolume within this target, and 7 organs at risk (OARs; other inserts) were contoured in the treatment planning system (TPS). Two adaptation sessions were performed in which the GTV was recontoured, first at 7 o'clock and then at 5 o'clock. The accumulated dose was exported from the TPS after delivery. Deformable vector fields were also exported to manually accumulate doses for comparison. For the pelvis phantom, synthetic Gaussian deformations were applied to the planning CT image to simulate organ changes. Two single-fraction adaptive plans were created based on the deformed planning CT and cone beam CT images acquired onboard the radiation therapy platform. A manual dose accumulation was performed after delivery using the exported deformable vector fields for comparison with the system-generated result. Results: All plans were successfully delivered, and the accumulated dose was both manually calculated and derived from the TPS. For the electron density phantom, the average mean dose differences in the GTV, boost volume, and OARs 1 to 7 were 0.0%, -0.2%, 92.0%, 78.4%, 1.8%, 1.9%, 0.0%, 0.0%, and 2.3%, respectively, between the manually summed and platform-accumulated doses. The gamma passing rates for the 3-dimensional dose comparison between the manually generated and TPS-provided dose accumulations were >99% for both phantoms. Conclusions: This study demonstrated agreement between manually obtained and TPS-generated accumulated doses in terms of both mean structure doses and local 3-dimensional dose distributions. Large disagreements were observed for OAR1 and OAR2 defined on the electron density phantom due to OARs having lower deformation priority over the target in addition to artificially large changes in position induced for these structures fraction-by-fraction. The tests applied in this study to a commercial platform provide a straightforward approach toward the development of routine quality assurance of dose accumulation in online adaptation.

First treatments for Lattice stereotactic body radiation therapy using magnetic resonance image guided radiation therapy.

Two abdominal patients were treated with Lattice stereotactic body radiation therapy (SBRT) using magnetic resonance guided radiation therapy (MRgRT). This is one of the first reported treatments of Lattice SBRT with the use of MRgRT. A description of the treatment approach and planning considerations were incorporated into this report. MRgRT Lattice SBRT delivered similar planning quality metrics to established dosimetric parameters for Lattice SBRT. Increased signal intensity were seen in the MRI treatments for one of the patients during the course of treatment.

DANSR: A Tool for the Detection of Annotated and Novel Small RNAs.

Existing small noncoding RNA analysis tools are optimized for processing short sequencing reads (17-35 nucleotides) to monitor microRNA expression. However, these strategies under-represent many biologically relevant classes of small noncoding RNAs in the 36-200 nucleotides length range (tRNAs, snoRNAs, etc.). To address this, we developed DANSR, a tool for the detection of annotated and novel small RNAs using sequencing reads with variable lengths (ranging from 17-200 nt). While DANSR is broadly applicable to any small RNA dataset, we applied it to a cohort of matched normal, primary, and distant metastatic colorectal cancer specimens to demonstrate its ability to quantify annotated small RNAs, discover novel genes, and calculate differential expression. DANSR is available as an open source tool.

LITE SABR M1: Planning design and dosimetric endpoints for a phase I trial of lattice SBRT.

PURPOSE: Lattice stereotactic body radiation therapy (SBRT) is a form of spatially fractionated radiation therapy (SFRT) using SBRT methods. This study reports clinical dosimetric endpoints achieved for Lattice SBRT plans delivering 20 Gy in 5 fractions to the periphery of a tumor with a simultaneous integrated boost (SIB) of 66.7 Gy, as part of a prospective Phase I clinical trial (NCT04133415). Additionally, it updates previously reported planning and delivery techniques based on extended experience with a broader patient population. METHODS: Patients were enrolled on a single-arm phase I trial conducted between November 2019 and August 2020. Eligibility was restricted to tumors >4.5 cm in the largest dimension. Characteristic SFRT dose gradients were achieved using a lattice of 1.5 cm diameter spheres spaced within the GTV in a regular pattern, with peak-to-valley dose varying from 66.7 Gy to 20 Gy within 1.5 cm. Organ-at-risk (OAR) sparing followed AAPM TG101 recommendations for 5-fraction SBRT. RESULTS: Twenty patients (22 plans) were enrolled on study, with one additional plan treated off study. All OAR and target coverage planning objectives were achieved, with the exception of a single small bronchus. Conformity of the 20 Gy isodose line significantly improved over the course of the study. The majority (85.2%) of treatment fractions were delivered in a 30 minutes timeslot, with 4 (3.5%) exceeding a total treatment time of 40 minutes. CONCLUSION: Lattice SBRT planning techniques produce consistent and efficient treatment plans. Refined techniques described here further improve the quality of the planning technique.

Technical Report: Development and Implementation of an Open Source Template Interpretation Class Library for Automated Treatment Planning.

PURPOSE: Widespread implementation of automated treatment planning in radiation therapy remains elusive owing to variability in clinic and physician preferences, making it difficult to ensure consistent plan parameters. We have developed an open-source class library with the aim to improve efficiency and consistency for automated treatment planning in radiation therapy. METHODS AND MATERIALS: An open-source class library has been developed that interprets clinical templates within a commercial treatment planning system into a treatment plan for automated planning. This code was leveraged for the automated planning of 39 patients and retrospectively compared with the 78 clinically approved manual plans. RESULTS: From the initial 39 patients, 74 of 78 plans were successfully generated without manual intervention. The target dose was more homogeneous for automated plans, with an average homogeneity index of 3.30 for manual plans versus 3.11 for automated plans (P = .107). The generalized equivalent uniform dose (gEUD) was decreased in the femurs and rectum for automated plans, with a mean gEUD of 3746 cGy versus 3338 cGy (P ≤ 0.001) and 5761 cGy versus 5634 cGy (P ≤ 0.001) for the femurs and rectum, respectively. Dose metrics for the bladder and rectum (V6500 cGy and V4000 cGy) showed recognizable but insignificant improvements. All automated plans delivered for quality assurance passed a gamma analysis (>95%), with an average composite pass rate of 99.3% for pelvis plans and 98.8% for prostate plans. Deliverability parameters such as total monitor units and aperture complexity indicated deliverable plans. CONCLUSIONS: Prostate cancer and pelvic node radiation therapy can be automated using volumetric modulated arc therapy planning and clinical templates based on a standardized clinical workflow. The class library developed in this study conveniently interfaced between the plan template and the treatment planning system to automatically generate high-quality plans on customizable templates.

An atlas-guided automatic planning approach for rectal cancer intensity-modulated radiotherapy.

We try to develop an atlas-guided automatic planning (AGAP) approach and evaluate its feasibility and performance in rectal cancer intensity-modulated radiotherapy. The developed AGAP approach consisted of four independent modules: patient atlas, similar patient retrieval, beam morphing (BM), and plan fine-tuning (PFT) modules. The atlas was setup using anatomy and plan data from Pinnacle auto-planning (P-auto) plans. Given a new patient, the retrieval function searched the top similar patient by a generic Fourier descriptor algorithm and retrieved its plan information. The BM function generated an initial plan for the new patient by morphing the beam aperture from the top similar patient plan. The beam aperture and calculated dose of the initial plan were used to guide the new plan optimization in the PFT function. The AGAP approach was tested on 96 patients by the leave-one-out validation and plan quality was compared with the P-auto plans. The AGAP and P-auto plans had no statistical difference for target coverage and dose homogeneity in terms ofV100%(p = 0.76) and homogeneity index (p = 0.073), respectively. The CI index showed they had a statistically significant difference. But the ΔCI was both 0.02 compared to the perfect CI index of 1. The AGAP approach reduced the bladder mean dose by 152.1 cGy (p < 0.05) andV50by 0.9% (p < 0.05), and slightly increased the left and right femoral head mean dose by 70.1 cGy (p < 0.05) and 69.7 cGy (p < 0.05), respectively. This work developed an efficient and automatic approach that could fully automate the IMRT planning process in rectal cancer radiotherapy. It reduced the plan quality dependence on the planner experience and maintained the comparable plan quality with P-auto plans.

Development of a storage phosphor imaging system for proton pencil beam spot profile determination.

PURPOSE: Accurate two-dimensional (2D) profile measurements at submillimeter precision are necessary for proton beam commissioning and periodic quality assurance (QA) purposes and are currently performed at our institution with a commercial scintillation detector (Lynx PT) with limited means for independent checks. The purpose of this work was to create an independent dosimetry system consisting of an in-house optical scanner and a BaFBrI:Eu2+ storage phosphor dosimeter by: (a) determining the optimal settings for the optical scanner, (b) measuring 2D proton spot profiles with the storage phosphors, and (c) comparing them to similar measurements using a commercial scintillation detector. METHODS: An in-house 2D laboratory optical scanner was constructed and spatially calibrated for accurate 2D photostimulated luminescence (PSL) dosimetry. Square 5 × 5 cm2 BaFBrI:Eu2+ dosimeter samples were uniformly irradiated with line scans performed to determine the physical and electronic scanner settings resulting in the highest signal-to-noise ratios (SNR) at a sub-millimeter spatial resolution. The resultant spatial resolution of the scanner was then quantitatively assessed by measuring (a) line pairs on a standard X-ray lead bar phantom and (b) modulation transfer functions. Following this, 2D proton spot profiles from a Mevion S250i Hyperscan proton unit were obtained at 1, 10, 20, 30, 40, and 50 monitor unit (MU) settings at maximum energy (E0  = 227.1 MeV) and compared to baseline profiles from a commercial scintillation detector, where 1 MU is calibrated to deliver 1 Gy absolute proton dose-to-water under reference conditions, that is, 41 × 41 proton spots uniformly spaced by 0.25 cm within a 10 × 10 cm2 square field size at maximum energy (227.1 MeV) in water at depth of 5 cm at isocenter. The dosimetric system's sensitivities to (a) ±1 mm positional shifts and (b) ±0.3 mm beam lateral spread changes were quantitatively evaluated through a Gaussian fitting of the crossline and inline plots of the respective artificially shifted beam profiles. RESULTS: The physical scanner settings of (a) Δτ = 27 ms time interval between data samples, (b) vx  = 1.235 cm/s scanning speed, (c) 1% laser transmission (0.02 mW power) and (d) (Δx, Δy) = (0.33, 0.50 mm) pixel sizes with electronic settings of (a) 300 microseconds time constant, (b) normal dynamic reserve, (c) 24 dB/oct low pass filter slope, and (d) 160 Hz chopping frequency resulted in the highest SNR while maintaining sub-millimeter spatial resolution. The BaFBr0.85 I0.15 :Eu2+ storage phosphor dosimeters were linear from 1 to 50 MU and their profiles did not saturate up to 150 MU. The scanner was able to detect lateral displacements of ±1 mm in both the crossline and inline directions and ±0.3 mm beam spread changes that were artificially introduced by varying the incident proton energy. Specific to our proton unit, proton energy changes of ±1 MeV can also be detected indirectly via beam spread measurements. CONCLUSION: Our combined dosimetric system including an in-house laboratory optical scanner and reusable BaFBr0.85 I0.15 :Eu2+ storage phosphors demonstrated a sufficient spatial resolution and dosimetric accuracy to support its use as an independent proton spot measurement dosimeter system. Its wide dynamic range allows for other versatile applications such as proton halo measurements.

Dynamics and extensional rheology of polymer-surfactant association complexes.

Understanding and characterizing the influence of polymers and surfactants on rheology, application, and processing is critical for designing complex fluid formulations for enhanced oil recovery, pharmaceuticals, cosmetics, foods, inks, agricultural sprays, and coatings. It is well-established that the addition of anionic surfactant like sodium dodecyl sulfate (SDS) to an aqueous solution of an oppositely-charged or uncharged polymer like poly(ethylene oxide) (PEO) can result in the formation of the polymer-surfactant association complexes (P0S-ACs) and a non-monotonic concentration-dependent variation in zero shear viscosity. However, the extensional rheology response of polymer-surfactant mixtures remains relatively poorly understood, partially due to characterization challenges that arise for low viscosity, low elasticity fluids, even though the response to strong extensional flows impacts drop formation and many processing operations. In this article, we use the recently developed dripping-onto-substrate (DoS) rheometry protocols to characterize the pinching dynamics and extensional rheology response of aqueous P0S- solutions formulated with PEO (P0) and SDS (S-), respectively. We find the PEO-SDS mixtures display a significantly weaker concentration-dependent variation in the extensional relaxation time, filament lifespan, and extensional viscosity values than anticipated by the measured shear viscosity.

Characterization of Dosimetric Differences in Strut-Adjusted Volume Implant Treatment Plans Calculated With TG-43 Formalism and a Model-Based Dose Calculation Algorithm.

PURPOSE: To comprehensively characterize dosimetric differences between calculations with a commercial model-based dose calculation algorithm (MBDCA) and the TG-43 formalism in application to accelerated partial breast irradiation (APBI) with the strut-adjusted volume implant (SAVI) applicator. METHODS: Dose for 100 patients treated with the SAVI applicator was recalculated with an MBDCA for comparison to dose calculated via TG-43. For every pair of dose calculations, dose-volume histogram (DVH) metrics including V90%, V95%, V100%, V150%, and V200% for the PTV_EVAL were compared. Features were defined for each case including (1) applicator size, (2) ratio between PTV_EVAL contour and 1-cm rind surrounding SAVI applicator, (3) ratio between dwell time in central catheter and total dwell time, and (4) mean computed tomography (CT) number within the lumpectomy cavity. Wilcoxon rank sum tests were performed to test whether treatment plans could be stratified according to feature values into groups with statistically significant dosimetry differences between MBDCA and TG-43. RESULTS: For all DVH metrics, differences between TG-43 and MBDCA calculations were statistically significant (P < .05). Minimum (maximum) relative percent differences between the MBDCA and TG-43 for V90%, V95%, and V100% were -2.1% (0.1%), -3.1% (-0.1%), and -5.0% (-0.5%), respectively. The median relative percent difference in mean PTV_EVAL dose between the MBDCA and TG-43 was -3.9%, with minimum (maximum) difference of -6.5% (-1.8%). For V90%, V95%, and V100%, plan quality worsened beyond defined thresholds in 26, 23, and 31 cases with no instances of coverage improvement. Features 1, 2, and 4 were shown to be able to stratify treatment plans into groups with statistically significant differences in dosimetry metrics between MBDCA and TG-43. CONCLUSIONS: Investigated dose metrics for SAVI treatments were found to be systematically lower with MBDCA calculation in comparison to TG-43. Plans could be stratified according to several features by the magnitude of dosimetric differences between these calculations.

Dual-storage phosphor proton therapy dosimetry: Simultaneous quantification of dose and linear energy transfer.

PURPOSE: To investigate the feasibility of using the high Zeff storage phosphor material BaFBrI:Eu2+ in conjunction with the low Zeff storage phosphor material KCl:Eu2+ for simultaneous proton dose and linear energy transfer (LET) measurements by (a) measuring the fundamental optical and dosimetric properties of BaFBrI:Eu2+ , (b) evaluating its compatibility in being readout simultaneously with KCl:Eu2+ dosimeters, and (c) modeling and validating its LET dependence under elevated proton LET irradiation. METHODS: A commercial BaFBrI:Eu2+ storage phosphor detector (Model ST-VI, Fujifilm) was characterized with energy dispersive x-ray spectroscopy (EDS) analysis to obtain its elemental composition. The dosimeters were irradiated using both a Mevion S250 proton therapy unit (at the center of a spread-out Bragg peak, SOBP) and a Varian Clinac iX linear accelerator with the latter being a low LET irradiation. The photostimulated luminescence (PSL) emission spectra, excitation spectra, and luminescent lifetimes of the detector were measured after proton and photon irradiations. Dosimetric properties including dose linearity, dose rate dependence, radiation hardness, temporal, and readout stabilities were studied using a laboratory optical reader after proton irradiations. In addition, its proton energy dependence was analytically modeled and experimentally validated by irradiating the detectors at various depths within the SOBP (Range: 15.0 g/cm2 , Modulation: 10.0 g/cm2 ). RESULTS: The active detector composition for the high Zeff storage phosphor detector was found to be BaFBr0.85 I0.15 :Eu2+ . The BaFBr0.85 I0.15 :Eu2+ material's excitation and emission spectra were in agreement under proton and photon irradiations, with peaks of 586 ± 1 nm and 400 ± 1 nm, respectively, with a full width at half maximum (FWHM) of 119 ± 3 nm and 30 ± 2 nm, respectively. As dosimeter response under photon irradiation is generally believed to be free from LET effect, these results suggest LET independence of charge storage center types resulted from ionizing radiations. There is sufficient spectral overlaps with KCl:Eu2+ dosimeters allowing both dosimeters to be readout under equivalent readout conditions, that is, 594 nm stimulation and 420 nm detection wavelengths. Its PSL characteristic lifetime was found to be less than 5 microseconds which would make it suitable for fast 2D readout post irradiation. Its 420 nm emission band intensity was found to be linear up to 10 Gy absolute proton dose under the same irradiation conditions, dose rate independent, stable in time and under multiple readouts, and with high radiation hardness under cumulative proton dose histories up to 200 Gy as tested in this study. BaFBr0.85 I0.15 :Eu2+ showed significant proton energy-dependent dose under-response in regions of high LET which could be modeled by stopping power ratio calculations with an accuracy of 3% in low LET regions and a distance-to-agreement (DTA) of 1 mm in high LET regions (>5 keV/µm). CONCLUSION: We have proven the feasibility of dual-storage phosphor proton dosimetry for simultaneous proton dose and LET measurements. BaFBr0.85 I0.15 :Eu2+ has shown equally excellent dosimetry performance as its low Zeff complement KCl:Eu2+ with distinctive LET dependence merely as a result of its higher Zeff . These promising results pave the way for future studies involving simultaneous proton dose and LET measurements using this novel approach.

Technical Note: Cumulative dose modeling for organ motion management in MRI-guided radiation therapy.

PURPOSE: To develop a method for continuous online dose accumulation during irradiation in MRI-guided radiation therapy (MRgRT) and to demonstrate its application in evaluating the impact of internal organ motion on cumulative dose. METHODS: An intensity-modulated radiation therapy (IMRT) treatment plan is partitioned into its unique apertures. Dose for each planned aperture is calculated using Monte Carlo dose simulation on each phase of a four-dimensional computed tomography (4D-CT) dataset. Deformable image registration is then performed both (a) between each frame of a cine-MRI acquisition obtained during treatment and a reference frame, and (b) between each volume of the 4D-CT phases and a reference phase. These registrations are used to associate each cine image with a 4D-CT phase. Additionally, for each 4D-CT phase, the deformation vector field (DVF) is used to warp the pre-calculated dose volumes per aperture onto the reference CT dataset. To estimate the dose volume delivered during each frame of the cine-MRI acquisition, we retrieve the pre-calculated warped dose volume for the delivered aperture on the associated 4D-CT phase and adjust it by a rigid translation to account for baseline drift and instances where motion on the cine image exceeds the amplitude observed between 4D-CT phases. RESULTS: The proposed dose accumulation method is retrospectively applied to a liver cancer case previously treated on an MRgRT platform. Cumulative dose estimated for free-breathing and respiration-gated delivery is compared against dose calculated on static anatomy. In this sample case, the target minimum dose and D 98 varied by as much as 5% and 7%, respectively. CONCLUSION: We demonstrate a technique suitable for continuous online dose accumulation during MRgRT. In contrast to other approaches, dose is pre-calculated per aperture and phase and then retrieved based on a mapping scheme between cine MRI and 4D-CT datasets, aiming at reducing the computational burden for potential real-time applications.

Direct tumor visual feedback during free breathing in 0.35T MRgRT.

To present a tumor motion control system during free breathing using direct tumor visual feedback to patients in 0.35 T magnetic resonance-guided radiotherapy (MRgRT). We present direct tumor visualization to patients by projecting real-time cine MR images on an MR-compatible display system inside a 0.35 T MRgRT bore. The direct tumor visualization included anatomical images with a target contour and an auto-segmented gating contour. In addition, a beam-status sign was added for patient guidance. The feasibility was investigated with a six-patient clinical evaluation of the system in terms of tumor motion range and beam-on time. Seven patients without visual guidance were used for comparison. Positions of the tumor and the auto-segmented gating contour from the cine MR images were used in probability analysis to evaluate tumor motion control. In addition, beam-on time was recorded to assess the efficacy of the visual feedback system. The direct tumor visualization system was developed and implemented in our clinic. The target contour extended 3 mm outside of the gating contour for 33.6 ± 24.9% of the time without visual guidance, and 37.2 ± 26.4% of the time with visual guidance. The average maximum motion outside of the gating contour was 14.4 ± 11.1 mm without and 13.0 ± 7.9 mm with visual guidance. Beam-on time as a percentage was 43.9 ± 15.3% without visual guidance, and 48.0 ± 21.2% with visual guidance, but was not significantly different (P = 0.34). We demonstrated the clinical feasibility and potential benefits of presenting direct tumor visual feedback to patients in MRgRT. The visual feedback allows patients to visualize and attempt to minimize tumor motion in free breathing. The proposed system and associated clinical workflow can be easily adapted for any type of MRgRT.

Quantitative proton radiation therapy dosimetry using the storage phosphor europium-doped potassium chloride.

PURPOSE: To (a) characterize the fundamental optical and dosimetric properties of the storage phosphor europium-doped potassium chloride for quantitative proton dosimetry, and (b) investigate if its dose radiation response can be described by an analytic radiation transport model. METHODS: Cylindrical KCl:Eu2+ dosimeters with dimensions of 6 mm diameter and 1 mm thickness were fabricated in-house. The dosimeters were irradiated using both a Mevion S250 passive scattering proton therapy system and a Varian Clinac iX linear accelerator. Photostimulated luminescence (PSL) emission spectra, excitation spectra, and luminescence lifetimes were measured for both proton and photon irradiations. Dosimetric properties including radiation hardness, dose linearity, signal stabilization, dose rate sensitivity, and energy dependence were studied using a laboratory optical reader after irradiations. The dosimeters were modeled using physical quantities including mass stopping powers in the storage phosphor and water for a given proton beam, and mass energy absorption coefficients and massing stopping powers in detector and water for a given photon beam. RESULTS: KCl:Eu2+ exhibited optical emission and stimulation peaks at 421 and 560 nm, respectively, for both proton and photon irradiations, enabling postirradiation readouts using a visible light source while detecting the PSL using a photomultiplier tube. KCl:Eu2+ showed a linear response from 0 to 8 Gy absorbed dose-to-water, a large dynamic range up to 60 Gy, dose-rate independence measured from 83 to 500 MU/min, and a PSL lifetime of <5 ms that is sufficiently short for supporting rapid scanning in a two-dimensional geometry. KCl:Eu2+ was highly reusable with only a slight signal decrease of ~3% at accumulated doses over 100 Gy, which could be managed by a periodic recalibration. The detected PSL signal strength of the dosimeter in the proton field had been calculated accurately to a maximum discrepancy of 2% using known physical quantities along with its prior signal strength as measured in a photon field at the same dose-to-water. This discrepancy might be attributed to an under-response due to linear energy transfer (LET) effect. However, comparisons of depth-dose measurements in a spread-out Bragg peak (SOBP) field with a parallel-plate ionization chamber showed no clear evidence of LET effects. Furthermore, range measurements agreed with ionization chamber measurements to within 1 mm. CONCLUSIONS: KCl:Eu2+ showed linear response over a large dynamic range for proton irradiations and reliably reproduced SOBP measurements as measured by ionization chambers. Its relatively low atomic number of 18 and near LET independence make it suited for quantitative proton dosimetry. In addition, its high radiation hardness means that it can be reused numerous times. Any potential measurement artifacts encountered in complex irradiation conditions should be able to be corrected for using known physical quantities.

Automated and robust beam data validation of a preconfigured ring gantry linear accelerator using a 1D tank with synchronized beam delivery and couch motions.

PURPOSE: To develop an efficient and automated methodology for beam data validation for a preconfigured ring gantry linear accelerator using scripting and a one-dimensional (1D) tank with automated couch motions. MATERIALS AND METHODS: Using an application programming interface, a program was developed to allow the user to choose a set of beam data to validate with measurement. Once selected the program generates a set of instructions for radiation delivery with synchronized couch motions for the linear accelerator in the form of an extensible markup language (XML) file to be delivered on the ring gantry linear accelerator. The user then delivers these beams while measuring with the 1D tank and data logging electrometer. The program also automatically calculates this set of beams on the measurement geometry within the treatment planning system (TPS) and extracts the corresponding calculated dosimetric data for comparison to measurement. Once completed the program then returns a comparison of the measurement to the predicted result from the TPS to the user and prints a report. In this work lateral, longitudinal, and diagonal profiles were taken for fields sizes of 6 × 6, 8 × 8, 10 × 10, 20 × 20, and 28 × 28 cm2 at depths of 1.3, 5, 10, 20, and 30 cm. Depth dose profiles were taken for all field sizes. RESULTS: Using this methodology, the TPS was validated to agree with measurement. All compared points yielded a gamma value less than 1 for a 1.5%/1.5 mm criteria (100% passing rate). Off axis profiles had >98.5% of data points producing a gamma value <1 with a 1%/1 mm criteria. All depth profiles produced 100% of data points with a gamma value <1 with a 1%/1 mm criteria. All data points measured were within 1.5% or 2 mm distance to agreement. CONCLUSIONS: This methodology allows for an increase in automation in the beam data validation process. Leveraging the application program interface allows the user to use a single system to create the measurement files, predict the result, and then compare to actual measurement increasing efficiency and reducing the chance for user input errors.